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ABSTRACT 

After studying the orientation of a warp generated by a companion satellite, we show 
that the Galactic Warp would be oriented in a different way if the Magellanic Clouds 
were its cause. We have treated the problem analytically and complemented it with 
numerical N-Body simulations. 
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1 INTRODUCTION 

The disk of the Milky Way is remarkably flat out to 10 
kpc, where it starts to bends in opposite directions in the 
southern and northern pa rts. T he cause of it is still a puzzle: 
for a review, see Binney (1992). 

One possibility is that the Magellanic Clouds distort the 
Galactic disk in the observed way. The fact that the direction 
of the maximum warping lies very close to the galactocentric 
longitude of the Clouds makes this hypothesis tempting. 

The problem with this scenario is that the Clouds are 
not massive enough to generate the warp amplitudes that we 
observe at their present dist ance. This w a s noticed b y the 
first researches in this field (Burke 1957; Kerr 1957), and 
later by Hunter & Toomre (1969). A remedy which might 
allow this scenario to work was to suppose that the Clouds 
are currently not at the pericenter of their orbit, so that 
they have been much closer to the Galactic disk in a pre- 
vious epoch (~ 20 kpc is what was needed). However, later 
work by Murai & Fujimoto (1986) determined the orbit of 
the Clouds, and proved that the Clouds are actually at their 
pericenter, with an apocenter close to 100 kpc, so the prob- 
lem of the small amplitude still remains if the Clouds are to 
be blamed as the cause of the Galactic warp. 

Recently a mechanism for amplify ing t he effect of a 
satellite has been proposed by Weinberg ( 1998 ) . He describes 
a calculation in which a disk galaxy surrounded by a dark 
halo is perturbed by a massive satellite, similar to the LMC. 
By means of a linear perturbation analysis, he follows the 
perturbation (wake) created by the satellite in the halo, in- 
cluding its self-gravity. He finds that the torque exerted on 
the disk is several times larger than that due directly to 
the satellite: the latter is amplified because (i) the satellite- 
induced wake in the halo itself exerts a torque, roughly in 
phase with that from the satellite; and (ii) the wake itself 



further perturbs the halo, resulting in a torque that is larger 
again. Under circumstances in which the satellite orbital fre- 
quency is close to the natural precession frequency of the 
disk ( i.e., the warp mode frequency of Sparke & Casertano 
(1988)), a sig nifican t amplitude can result. A calculation by 



Lynden-Bell (1985) of a similar scenario gives comp arabl e 



results, as does a simple model described by Kuijken (1997) 



In this paper, we focus on the orientation of a warp 
generated by a massive orbiting satellite with less emphasis 
on the amplitude of the warp. In §2 we discuss a simple 
analytic model in which the disk and halo are rigid: this 
establishes the baseline response of a disk to satellite tides, 
and its orientation with respect to the satellite orbital plane. 
As we show, this orientation is different from that of the 
Galactic warp to the LMC orbit. §3 contains a description 
of the N-body code used, and §4 to §6 the results of our 
N-body simulations, showing that the orientation problem 
remains. In §7 we give our conclusions. 



2 ANALYTIC RESULTS WITH A SIMPLIFIED 
MODEL 

A simple model serves as a reference for the warp response 
of the disk to an orbiting satellite. Consider a rigid disk, 
embedded in a rigid halo potential, and subjected to the 
potential of an orbiting satellite. The evolution of the disk is 
governed by the combined torque from halo and satellite. A 
stellar or gaseous disk is floppy, and so will warp when tilted, 
since it is not able to generate the stresses that would be 
required to keep it flat; however the overall re-alignment of 
the disk angular momentum should be comparable between 
the rigid and floppy cases. 

The angles used in this paper related to the satellite, 
and the definition of our coordinate system are illustrated 
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Figure 1. Definition of the coordinate system, satellite angles, 
and orientation of the disk of the Galaxy. The disk lies on the 
z = plane. The sun is on the left, and the galactic plane is 
viewed from the South Galactic Pole so that the disk rotates 
counterclockwise (indicated by an arrow). 

in Figure |l|. The tilting of the disk is measured by the angle 
between the z axis and the angular momentum of the disk. 
The longitude of this vector is the same as the longitude of 
the maximum of the warp when looking at the disk from the 
North Galactic Pole. 

The Lagrangian for a rigidly spinning, axisymmetric ob- 
ject is 

C = \h{6 2 + <j> 2 sin 2 e) + ±I a {<j>cos9 + i>) 2 - V{9, <j>) (1) 

where (9, <f), ip) are the Euler angles, and 73 and h are the 
moments of inertia of the object about its symmetry axis and 
about orthogonal directions. V is the potential energy of the 
body in the halo plus satellite potential. The i^-equation of 
motion leads to the conserved quantity S = l3(4>cos9 + ip), 
the spin, and the other two equations of motion then become 

■o ■ dV 

h6- h4> 2 sin cos + S<j>sinO + =0 (2) 



i)0 



and 
d 



h 



dV 



(3) 



For small deviations from the equator (8 — 0), we can ex- 
pand these equations in terms of x = sin 9 cos cj> ~ 6 cos <f>, 
y = sin#sin</> ~ #sin0. In these terms the equations of 
motion become 
dV 
dx 
dV 
dy 

For small x, y, the potential energy of the disk due to the 



h'x + Sy + 
hy — Sx 



0, 



0. 



(4) 
(5) 



i y ), and that 
will be — Vsfsin 2 9s — 



flattened halo will have the form 
due to the satellite at position 9s,<j)s 
x sin 29s cos 0s — y sin 29s sin (j>s) where Vh and Vs are con- 
stants representing the strengths of the halo torque and of 
the quadrupole of the tidal field from the satellite, respec- 
tively. Hence we find 



h'x + Sy + Vhx + Vs sin 29s cos 4>s = 0, 
hy — Sx + VhV + Vs sin 29s sin 4>s = 0. 



(6) 
(7) 



If furthermore the satellite orbit is circular and polar in the 
x — plane, 9s = £lst, 0s — 90, and the solution to the 
equations of motion is 



2Q S S 

x — — — — Vscos2f2st; 



V = 



4fin§ 



V u 



Vssin 2Q s t 



(8) 



plus free precession and nutation terms, where A = (Vh — 
4Jir2l) 2 - 40,1s 2 . (A more general quasiperiodic satellite 
orbit yields a solution which can be written as a sum of 
such terms.) Notice that the satellite provokes an elliptical 
precession about the halo symmetry axis, with axis ratio 
dependent on the halo flattening and on the satellite orbit 
frequency. For example, for an exponential disk of mass M, 
scale length h and with a flat rotation curve of amplitude 
v, I 3 = 2h = 6Mh 2 and S = 2hvM. For such a disk in 
a spherical (or absent) halo (Vh = 0), a satellite orbiting 
at radius rs has frequency fis = v/r, and hence the axis 
ratio of the forced precession is (x : y) = rs/Sh. Hence the 
response of the disk to a distant satellite is mainly to nod 
perpendicular to the satellite orbit plane. This result can be 
understood as the classic orthogonal response of a gyroscope 
to an external torque: a distant satellite has a sufficiently low 
orbital frequency that the disk responds as if the torque were 
static. 

For a slightly flattened potential of the form ^v 2 ln[7? 2 -|- 
(2/(1 — e)) 2 ], Vh = Mv 2 e. With non-zero e, the axis ratio of 
the precession cone becomes [(4/i/r S )/(e - 12/i 2 /t|)]: again 
the oscillation in x is larger than that in y except for very 
flattened halos. 

The amplitudes generated by tidal perturbation from 
a satellite such as the LMC are small, less than a degree. 
The largest amplitude of oscillation is in the y-direction. 
The potential energy of the disk due to the tidal field of the 
satellite can be shown to be (see Appendix) 



Vs = 



3GM s h 
24 



Hence equation 
amplitude of 

9 GM S h 
8 v 2 rs rs 



0.09 



(9) 

yields, to leading order in h/rs, an x- 



(10) 



for the LMC (orbital radius about 50kpc, and rs/h ~ 15). 
This number increases only slightly (a factor 2) for halo 
flattenings up to 0.2 (see Figure Eh. 

It is clear from this calculation that simple tidal tilting 
of a disk by an LMC-like satellite does not provide a good 
model for the warping in the Galaxy, because the orientation 
of the warp is not perpendicular to the orbital plane of the 
LMC. This constraint is independent of the strength of the 
perturbation Vs- 

The amplitudes are much too small, but we have only 
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Figure 2. The oscillation of the axis of a rigid exponential disk 
subjected to the tidal field of an orbiting satellite. The amplitude 
is calculated assuming a satellite of mass 1.5 X 1O 1O M , orbiting 
at radius 54kpc in the z = plane. The direction of the tilt of 
the Galactic disk with respect to the Magellanic Clouds' orbital 
plane is indicated by the arrow. The dots mark the expected 
position of the disk axis given the current phase of the LMC 
orbit for (bottom to top) halo potential ellipticities e = (solid 
symbol), 0.05, 0.1, 0.15 and 0.2 (open circles). 



considered thy Lilling of a, rigid dink, and tire sitnattan ran 
change ' when Che floppiness of the disk is considered. 



3 SIMULATION DETAILS 

To test this scenario, and in particular to get beyond the 
rigid tilting considered above, we have performed some N- 
body simulations. We assume the halo to be a background 
potential which does not respond to the disk or the satellite. 

The N-body code used for this work models the disk as 
a set of concentric spinning rings embedded in a spherical, 
rigid halo. This description allows warps to be described, 
but not in-plane distortions of the disk such as bars or lop- 
sidedness. 



3.1 Initial conditions 

We have performed simulations with two types of disks: a 
rigid disk and a exponential disk. The rigid disk run tells us 
how good the analytic predictions are, and the exponential 
disk is used later for a more realistic approach. 

We use a King model for the halo, in order to obtain 
a reasonable flat rotation curve (Figure ^). This is accom- 
plished with a model of vPo/cto 2 = —6, a tidal radius of 44 
(200 kpc. for a 4.5 kpc scale-length disk), and mass of 10 
disk masses. 

Each of the disks is made of 1000 rings, each of them 
consisting of 36 particles. Various runs where made with 
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Figure 3. Rotation curve showing the contribution from the disk 
(dotted) and halo (dashed) to the total (solid) 



more rings and more particles per ring, without important 
changes in the results described below. 

The satellite is modelled as having a Plummer distribu- 
tion. To avoid relative movements of the galaxy with respect 
to the satellite we have used two satellites instead of one, 
symmetrically placed with respect to the centre of the halo- 
disk system. This causes the dipole term of the tidal field 
to be zero, avoiding relative movements of the galaxy with 
respect to the satellites. It is equivalent to only keeping the 
even-/ terms in the potential of the satellite, neglecting the 
dipole, I — 1, terms in the potential, and concentrating on 
the warp (which result from the quadrupole, I = 2 terms). 

The first run was made with a satellite in a circular 
orbit, to try to reproduce the predictions in §2. Later a non- 
circular orbit is considered, and the difference between both 
simulations analysed. The non-circular orbit has a pericen- 
ter at 50 kpc and apocenter at 100 kpc, consi stent with re- 
cent determin ation of the orbit of the Clouds (Lin, Jones & 



Klemola 1995 ). In the non-circular simulations the satellite 
starts at its apocenter at the beginning of the simulation, 
where the perturbation on the disk is the smallest possible. 

The units of the model translate to the Galaxy (disk 
scale-length 4.5 kpc, and the rotation velocity at 8.5 kpc of 
220 kms -1 ) as follows: length unit = 4.5 kpc, velocity unit 
= 340 kms -1 , time unit = 1.30 x 10 7 years, mass unit= 
1.20 x 1O 11 M0. With these numbers, the disk mass of our 
model is 6.1 x 10 10 Mq, and the satellite (LMC) has a mass 
of 1.5 x 10 10 Mpj, the biggest current mass estimate for the 
Clouds dSchommer et al. 1992 ) . In the coordinate system of 
the simulations, the z — is the disk plane, and the orbit of 
the satellite lies in the x — plane. 



3.2 Code used to evolve the system 

The disk is modelled as a system of pivoted spinning rings, 
fixed at the centre of the halo. Each ring is realized as 36 
azimuthally-spaced particles, and the p otential generated b y 
the rings is calculated with a tree code (Barnes & Hut 1986). 



The forces on the individual "ring-particles" are used to cal- 
culate the torque on each ring. The force exerted by a satel- 
lite on the ring particles was evaluated directly using the 
Plummer law. 

The Euler equations for rings and axisymmetric bodies 
can be rewritten in a form so that the time derivatives of 
the instantaneous angular velocities about the body axes are 
linear combinations of the angular velocities, torques and 
body normal vector components. This allows the derivation 
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Precession path for Rigid Disk 




Figure 4. Precession path followed by the rigid disk: x = 8 cos <j>, 
y = Bs'mifi. The dots indicate the disk's state when the satellite 
has the LMC's orbital phase. 



of a second order explicit time-centred leapfrog integration 
scheme that can be used to solve the coupled equations for 
t he rings and m ake it easy to merge with an N-body code 
(pubinski 200fj|). 



4 RIGID DISK 

As a first approach, we have evolved a rigid disk and anal- 
ysed its evolution under the influence of an orbiting satellite. 
The result of our simulation is in good agreement with the 
analytic predictions. The disk wobbles under the influence 
of the satellite, describing an ellipse elongated in the direc- 
tion perpendicular to the satellite's orbital plane. The path 
followed by the disk is plotted in Figure ^| where it can be 
seen that most of the time the maximum of the warp is lo- 
cated in the direction perpendicular to the orbital plane of 
the satellite (Z = 0° and I = 180°). The ellipse isn't as regu- 
lar as in Figure ^| for two reasons: the assumption that the 
disk is much smaller than the orbital radius of the satellite is 
not completely fulfilled; and there are some transient terms 
present because of the initial conditions of the simulation. 
This are also the cause for the precession ellipse of the disk 
not to be centred in the origin. 

The position of the warp when the satellite is at the 
location of the LMC is indicated by the dots in Figure ^, and 
the location of them resembles the predicted one in Figure ^ 
(for thaio = 0) remarkably well. 



5 EXPONENTIAL SELF-GRAVITATING DISK 

We now consider a more realistic disk: an exponential disk 
model, in which we have considered also the disk's self- 
gravity. The first thing that draws our attention in this 
simulation is a peak we see in the inclination at around 
6.5 scale-lengths. Simulations done with a different rotation 
curve showed that this peak occurs at the locations on the 
disk that satisfy Q s /w z = 2, 3, that are caused by reso- 
nances with the satellite's orbital frequency. This happens 
because the non-linear behaviour of the outer parts of the 
disk, where the assumption r s 3> rdisk is worse than it is 
further in. 

This is not the kind of warp we are looking for, due 
to the fact that it is the result of a satellite with a single 
frequency, and in the real case the eccentric orbit of the 
satellite will wash out this peak. Looking at the evolution 
of the disk it is clear that the warp looses its coherence at 
a radius about 4.5 scale-lengths (at larger radius the line 



of nodes winds up), so we will measure the warp properties 
considering that the disk finishes there. 

In the case of a floppy disk it is not straight-forward 
to define a single inclination and position angle. We have 
separated the disk into two components: the inner disk and 
the outer (warped) disk. The inner disk consists on the first 
2 scale- lengths, and remains practically flat along the simu- 
lation. The warping angle is then calculated as the angle 
between the inner and outer disk vectors. We have chosen 
to use the disk vectors and not the angular momentum, for 
example, not to penalise the outer less massive rings. The 
results presented here do not change significantly when the 
definition of the inner disk is altered. 

It has to be kept in mind that the warping angles quoted 
here are different than the maximum amplitude of the warp, 
who usually are larger by a factor not greater than 5. 

Using this method we obtain a plot similar to Figure [| 
for the exponential disk, which is shown in Figure ^. Only 
the path after t=f60 is shown, that is the moment when the 
disk behaviour reaches an equilibrium. 

Note that the predictions for the Galactic Warp don't 
really change with the floppiness of the disk: it is clearly close 
to I = 0°, as chapter §2 predicted, and not at I =~ 90°, as 
we observe it in the Galaxy. 



6 NON-CIRCULAR ORBIT, AND FLATTENED 
HALOS 

We also considered non-circular orbits, to allow for the fact 
that the orbit derived for the Clo uds has a pericenter of 5 
kpc and an apocenter of fOO kpc (Murai & Fujimoto f980). 



The changing radius of the satellite causes a fluctuating tidal 
field amplitude, which could be important for the dynamics 
of the disk. Here we show that in fact the effect does not 
change our conclusions materially. 

First, to have an idea of what to expect, we integrated 
the analytic equations of section §2 with a satellite in this 
kind of orbit. The result was, as before, that the disk's 
precession path was contained within an ellipse, elongated 
along the direction perpendicular to the satellite's plane. 
This causes the warp maxima to be most of the time close 
to the direction perpendicular to the satellite's orbit. 

We then performed simulations with this type of or- 
bits. The first thing we observe in these simulations is that 
the resonance peak we found in the circular orbit simulation 
has disappeared. Now the satellite doesn't have a single fre- 
quency, so the result is not surprising. The energy of the 
resonance gets distributed along different parts of the disk 
now, and no coherent pattern can be maintained across the 
disk, winding up the outer parts of the disk. When we look at 
the inner 4.5 scale-lengths as before, the precession pattern 
remains similar to the simulation with the circular orbit, so 
does the prediction of the warp's longitude at LMC's actual 
orbital phase. So our conclusions are not modified by the 
non-circularity of the orbit. 

The halos considered in all these simulations are spher- 
ical, which means that they don't contribute to the gener- 
ation of torques on the disk. We know that halos are not 
spherical, which creates a preferred plane in which the disk 
settles. Ellipticities of the order of 0.05 in the potential make 
the precession paths described before yet more elongated, 
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Figure 5. Warp orientation followed by the Exponential Disk: x = 9cos<f>, y = 8s'm<f>. The dots indicate the disk's state when the 
satellite has the LMC's orbital phase. 



which would make the chances of finding the warp maxima 
in the satellites' direction even more unlikely. 



7 CONCLUSION 

We show by means of analytic work and N-body simulations 
that the precession path of a warp generate by an orbiting 
satellite galaxy is elongated along the direction perpendic- 
ular to the satellite's orbital plane. Applying our result to 
the Milky Way, if the Galactic Warp is generated by the 
Magellanic Clouds, the direction of maximum amplitude of 
the warp would lie close to I ~ 0°, as compared to the ob- 
served direction of ~ 90° . Even if the halo's self-gravitating 
tidal response to the satellite amplifies the effect of the satel- 
lite (Weinberg 1998), this response will be mostly in phase 
with the satellite, and the alignment problems will persist. 
Possibly the Sagittarius dwarf galaxy, whose orbit lies at 
90° to that of the LMC, can b e the cause of the warp in- 
stead ([bata & Razoumov 1998). 



A limitation of the present work is that the halo has not 
been conside red as a live, self-gravitating component. It has 



been shown (Dubinski & Kuijken 1995; Nelson & Tremainc 



1995) that the back-reaction of the halo on a re-aligning 



disk can have important consequences. Such effects will be 
the subject of a further paper. 



APPENDIX A: POTENTIAL OF 
AXISYMMETRIC DISK DUE TO A SATELLITE 

The potential energy of an disk of surface density £(r) and 
in the gravitational field due to a satellite at position rs is 
given by 

V = - f d 2 rGS(r) Ms , . (Al) 
J \r-rs\ 
Choosing spherical coordinates for the satellite's posi- 
tion (see Figure |l|), and Cartesian coordinates in the disk 
plane so that the satellite has x = 0, we have 



V = -GM S 



T,dx dy (r-g — 2j/rs sin 9s + x 2 + y 2 ) 



-1/2 



(A2) 



Assuming that the disk is small compared to rs, we can 
expand the integrand in x and y. For an axisymmetric disk 
the second-order terms are the first ones that generate a 
potential gradient: they are 



V = 



GM S 



J T,dxdy [-i(l-3sin 2 <9 s )j/ 2 - \x 2 ] 



which results in 



V 



3GM s h . 2r . 
= — sin as + constant. 



2r§ 



(A3) 



(A4) 
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